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We consider relativistic, stationary, axisymmetric, polytropic, unconfined, 
perfect MHD winds, assuming their five lagrangian first integrals to be known. 
The asymptotic structure consists of field-regions bordered by boundary layers 
along the polar axis and at null surfaces, such as the equatorial plane, which have 
the structure of charged column or sheet pinches supported by plasma or mag- 
netic poloidal pressure. In each field-region cell, the proper current (defined here 



X 



as the ratio of the asymptotic poloidal current to the asymptotic Lorentz fac- 
tor) remains constant. Our solution is given in the form of matched asymptotic 
solutions separately valid outside and inside the boundary layers. An Hamilton- 
Jacobi equation, or equivalently a Grad-Shafranov equation, gives the asymptotic 
structure in the field-regions of winds which carry Poynting flux to infinity. An 
important consistency relation is found to exist between axial pressure, axial 
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current and asymptotic Lorentz factor. We similarly derive WKB-type analytic 
solutions for winds which are kinetic-energy dominated at infinity and whose 
magnetic surfaces focus to paraboloids. The density on the axis in the polar 
boundary column is shown to slowly fall off as a negative power of the logarithm 
of the distance to the wind source. The geometry of magnetic surfaces in all 
parts of the asymptotic domain, including boundary layers, is explicitly deduced 
in terms of the first-integrals. 



1. Introduction 

Any stationary axisymmetric non-relativistic, rotating, magnetized wind will collimate 
at large distances from the source, under perfect MHD conditions and polytropic thermody- 
namics (Heyvaerts & Norman 1989). Chiueh et al. (1991) showed that these results hold also 
for relativistic winds. We have recently extended our general analysis (Heyvaerts & Norman 
2002a) by presenting global analytic asymptotic solutions for non-relativistic winds, valid 
from the pole to the equator, assuming given first integrals. This paper extends our general 
analysis to relativistic winds. 

Flows which bring no Poynting flux to infinity are called kinetic winds. Their magnetic 
surfaces asymptote to paraboloids (Heyvaerts & Norman (1989), Chiueh et al. (1991)). If 
Poynting flux reaches infinity, the flow is called a Poynting jet. The magnetic surfaces then 
asymptotically approach cylinders or conical surfaces in which cylindrical ones are nested. 
We refer to a magnetic surface as being asymptotically conical if both r and z approach 
infinity on this magnetic surface, such that (z/r) approaches a finite value. 

We find that the asymptotic structure of relativistic winds consists of field-regions where 
the Lorentz force vanishes in the direction normal to magnetic surfaces. These regions are 
bounded by null surfaces where the magnetic field vanishes. In the vicinity of null surfaces, 
the plasma pressure, or possibly the poloidal magnetic pressure, is significant. The vicinity 
of the polar axis is also a boundary layer region. 

Relativistic winds are likely to be present near pulsars, active galactic nuclei, micro- 
quasars and gamma ray burst sources. Their structure has been first discussed assuming 
conical magnetic surfaces (Michel (1969), Goldreich & Julian (1970)). The cross field force 
balance has been considered by Okamoto (1975) and Heinemann & Olbert (1978). It has 
then been recognized that both non-relativistic and relativistic winds focus asymptotically 
(Blandford k Payne (1982), Heyvaerts & Norman (1989), Chiueh et al. (1991)) and that 
MHD forces may contribute to the wind acceleration. It is possible however that a confine- 



- 3- 



ment mechanism other than the action of the hoop stress should be operating (Spruit et al. 
(1997), Lugek & Bell (1997)). The general-relativistic dynamics of MHD winds for given field 
geometries has been discussed by Camenzind (1986a,b, 1989) and by Takahashi et al. (1990). 
Much work has also been devoted to the determination of the shape of the magnetic surfaces 
and to the wind dynamics. These studies considered the special-relativistic (Ardavan (1979), 
Camenzind (1987), Li (1993a,b)) as well as the general-relativistic transfield equation that 
describes force balance accross field lines (Mobarry & Lovelace (1986), Camenzind (1987), 
Nitta et al. (1991), Beskin & Par'ev (1993) Tomimatsu (1994), Koide et al. (2000)). 

Analytical models, usually involving some form of self-similarity, have been presented 
(Contopoulos (1994), Contopoulos (1995), Contopoulos & Lovelace (1994), Lovelace et al. 
(1991)). Li et al. (1992) extended the self-similar analysis of Blandford & Payne (1982) 
to relativistic winds. Their approach, as is the case here, does not rely on the force-free 
approximation. Their self-similarity requirement restricts the rotation profile of the wind 
source to be proportional to r -1 , imposes the wind source to be of an infinite extent and 
the return current flow to be an axial singularity. This limits the type of asymptotics 
for this model to circum-polar current carrying structures. The parabolic shape obtained 
by Li et al. (1992) is made possible by the fact that their wind source subtends infinite 
flux (Heyvaerts & Norman (2002a), appendix C). By contrast, their approach deals with 
the criticality and Alfven regularity conditions exactly, whereas in our approach it is only 
assumed that the set of first integrals is consistent with these conditions. Because the return 
current distribution is most important for determining the degree of collimation, we regard 
a non-self-similar approach as preferable for analyzing the asymptotics of the wind. We also 
have the possibility of analyzing a point source of wind which subtends finite total flux. 

Solutions for relativistic force-free flows have beeen obtained in cylindrical symmetry, 
either for given current profiles (Appl and Camenzind (1993a), Appl and Camenzind (1993b)) 
or for self-consistent plasma flows both confined (Beskin & Malyshkin (2000)) or unconfined 
(Nitta (1997)). Beskin et al. (1998) construct a solution by expanding in terms of the 
inverse of Michel's magnetization parameter. Cold relativistic winds have been analyzed 
by Bogovalov (1999) for the oblique split-monopole rotator. The question of the structure 
and formation of jets and their connection to disks have been dealt with numerically in 
many papers. Axisymmetric stationary solutions, enforcing regularity at the light cylinder, 
have been obtained by Contopoulos et al. (1999). Axisymmetric two-dimensional force-free 
stationary flows emitted by disks have been calculated by Fendt et al. (1993), Fendt & 
Camenzind (1996), Tsinganos & Bogovalov (1999) and in three dimensions by Krasnopolsky 
et al. (1999), Koide et al. (1999) and Nishikawa et al. (1997). The confinement of a star's 
flow by the wind from a disk has been studied by Bogovalov & Tsinganos (2001). Van 
Putten (1997) studied time-dependent simulations of magnetized relativistic jets. Bogovalov 
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(2001) finds that in relativistic winds the collimated region may reduce to a very small region 
trapping but little flux. Lovelace and collaborators considered the jet-disk interaction and 
have used the force free approximation to study the dynamics and focusing of relativistic 
jets (Lovelace (1976), Lovelace et al. (1987), Lovelace et al. (1991), Lovelace et al. (1993), 
Ustyugova et al. (2000)). The self-consistent response of the disk emitting the wind has also 
been considered in two dimensions (Bell & Lugek (1995), Kudoh et al. (1998), Koide et al. 

(2002) ), and in three dimensions (Matsumoto et al. 1996)). Non-stationary behaviour and 
an avalanche type of accretion flow results. 

Special analytical solutions valid in more or less extended regions of space have been 
obtained. Li (1993a) and Begelman & Li (1994) give some elements of our general solution 
below, although without the complete matching. In studying the structure of relativistic 
winds far from the light cylinder, Tomimatsu (1994) adopted a procedure similar in spirit 
to our present work. Tomimatsu & Takahashi (2003), in an elegant asymptotic analysis, find 
similar slow logarithmic wind and jet acceleration. Nitta (1995) developed a particular 
solution in the limit of winds with a very large mass flux to magnetic flux ratio and rigid 
rotation. He matched a solution valid in a cylindrical polar region, in which a strong current 
flows, with a conical solution in the nearby region. The force balance in the axial region 
is between the centrifugal force and the hoop stress, implying that this region is not much 
broader than the light cylinder. Some of our general asymptotic solutions also have such a 
mixed structure, although the force balance near the axis is different. Okamoto (1999) insists 
on the importance of obtaining solutions consistent with current closure and has shown that 
this implies that, in regions where the electric current returns to the wind source, magnetic 
lines should bend toward the equator instead of bending to the axis (Beskin & Okamoto 
2000; Okamoto 2000). The solutions presented below, although their poloidal lines curve to 
the polar axis in most of the plasma volume, comply with this requirement. 

It is the aim of this paper to provide a general analytical asymptotic solution for special- 
relativistic jets, assuming the five first integrals of the motion to be known. It is organized as 
follows. Section (2) recalls the basics of special-relativistic, stationary axisymmetric rotating 
MHD winds. Section (3) deals with the field-regions. The asymptotic form of the transfield 
equation in field-regions is established and reduced to a simple Hamilton-Jacobi equation 
which we solve in section (4). The solution in the polar boundary layer is then obtained in 
section (5), assuming that it encompasses little flux. Our solution applies both to Poynting 
jets and to kinetic winds by means of a WKB approximation. We match the field-solution 
to that which applies to the polar boundary layer. This specifies how the proper current 
around the polar axis varies with distance to the wind source. The case of a polar boundary 
layer supported by the poloidal magnetic pressure is also considered for completeness. In 
section (6) we similarly obtain and match to the field-region a solution valid in the vicinity 
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of a null magnetic surface, namely the equatorial plane of a magnetic structure with dipolar 
type of symmetry. In section (7) the shape of the magnetic surfaces is calculated, both 
in field-regions and in boundary layers, for both cases of asymptotic regime. Conclusions 
regarding the general properties of relativistic rotating MHD winds are presented in section 
(8). 



2. Axisymmetric Stationary Relativistic MHD Flows 

2.1. Notation and Definitions 

We now review relativistic MHD winds and establish our notations. We use cylindrical 
coordinates (r, 9, z). Unit vectors of the associated local frame of reference are e r , eg, e z . 
The notation R denotes the spherical radius. The unit normal vector to the poloidal field 
lines, oriented towards the polar axis, is ft and the unit tangent vector to them, oriented 
towards increasing z, is t. This vector makes an angle ip with its projection on the equatorial 
plane. Any vector field can be split into a poloidal (subscript P) and a toroidal (subscript 
6) part. Due to axisymmetry, the poloidal magnetic field, B P , can be expressed in terms of 
a magnetic flux function, a(r, z), such that 

- 1 da _ 1 da _ 

B P = — Tre r + -—e z . (1) 
r oz r or 

A magnetic surface is generated by rotating a field line about the polar axis. It is a surface 
of constant value of a(r,z). The magnetic flux through it is 2na. The flow, described here 
in the framework of special relativity, has a local Lorentz factor, 7 defined by: 

We denote by p the proper rest mass density, measured in the rest frame of the fluid. The 
momentum per unit proper mass is: 

u — ^fv (3) 



2.2. First Integrals 

A polytropic law is assumed, by which the proper gas pressure is related to the proper 
density by 

P = Q(a)p r (4) 
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where Q is constant following the fluid motion and T is the polytropic index. This relation 
may represent adiabatic or more complex thermodynamics. We define the function 



1 + r— i— ■ (5) 



which is also equal to (1 + f dP/pc 2 ) calculated at constant polytropic entropy Q. Denoting 

— * 

by p e the electric charge density, by j the electric current density and by M* the mass of 
the central object, the special-relativistic equation of motion can be written as (Goldreich 
& Julian 1970; Li 1993a): 

ipiv- V)( 7 ^) = —VP + jxB + p e E + 7P V (6) 

The relativistic form of the laws of mass conservation, isorotation, angular momentum conser- 
vation and Bernoulli are obtained as in the non-relativistic case and involve surface functions 
E, a, L, Q and Q. The equations which express these four laws are: 

-fpvp = a(a)B P , (7) 

^ve - rVt{a)) = a(a)B e /p, (8) 

2 GM, r£l(a)Be 

7e(c " " = m (10) 

Note that E includes the rest mass energy. The polytropic factor Q is a surface function 
because, by stationarity, flow surfaces are also magnetic surfaces. The rotation rate of the 
magnetic field, Q(a), which appears in equations (8) and (10) is defined in terms of the 
electric field by: 

E = -VL(a)Va (11) 

When Bp is directed away from the wind source, a is positive. Since the sense of the 
magnetic field is immaterial, it can be assumed that this is so at the positive polar axis. For 
positive a, a increases from pole to equator. We assume Q to be always positive. 



2.3. Bernoulli Equation 

The toroidal variables may be eliminated by using Eqs.(8) and (9). This gives: 

rB e = p ap — (12) 

/i a 2 4 - p 
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7^0 = - + — c 13 

We denote by / the quantity 

'—ir (14) 

Po 

The minus sign in Eq.(14) has been included to make / positive when a and Vl are. The 
physical total poloidal current is J P = —2-kI. Nevertheless we shall conveniently refer to / 
as the poloidal current. Since 7 depends on vg, the elimination of the toroidal variables in 
Eqs.(12) and (13) is not yet complete. These expressions can however be substituted in the 
Bernoulli equation (10) which can then be solved to obtain an expression of 7^ in terms of 
poloidal variables. This eventually gives the toroidal variables in terms of poloidal ones as: 

L(c 2 - GM*/R) - r 2 tlE 
rBd ~ ^° a (c 2 -GMjR)(l-^/p)-r^ (15) 

rv e = r^(l r^E-L(c 2 -GMjR) \ 

rVe riZ V + p r 2 nE (l-p a 2 ap)-LrW) 1 &j 

For smooth continuuous solutions, these expressions must be regular where their denominator 

vanishes, which implies that when 

c 2 - GMJR 

P = ^ a ^ c 2 -GMjR-n^ (17) 
the position (r, z) must be such that 



r2= L (c 2 -GM*/V^+^) 

n e 1 ' 

Another way to express the special density which appears in (17) is to insert in it the expres- 
sion (18) for the corresponding radius, so obtaining a relation between the value assumed by 
p at this special point and the first integrals: 

Poa 2 -Att^ (19) 



i{ P ) — E-m 

A complete elimination of toroidal variables from Eq.(10) can be achieved as follows. A first 
expression for 7 is obtained by substituting Eq.(12) in Eq.(10). An independent relation for 
7 results from its definition in Eq.(2), using Eqs.(7) and (13). Eliminating 7 between these 
two relations, an equation for p, or any other poloidal variable, is obtained. For given values 
of the first integrals, this relation, the relativistic Bernoulli equation, can be used to find 
this poloidal variable as a function of position along the magnetic surface a. Let us denote 
by D the variable: 

D = L \ ( 20 ) 
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and by g the function: 

9 = GM ' (21) 
c 2 yr 2 + z 2 

Simple algebra shows that the relativistic Bernoulli equation can be then written as: 

W -,)- DW -,)-w))' (e +1 ^) = 

(£ _ _ m) f + / _ LQ )I__ _ LQ (i _ 5 ) \ (22) 

This equation is satisfied at the relativistic Alfven point, as can be shown from Eqs.(18)- 
(19). Again, a smooth solution of Eq.(22) requires that the first integrals satisfy regularity 
conditions by passing critical points. 



2.4. Transfield Equation 



The transfield equation is the projection of the equation of motion on the normal to 
magnetic surfaces. It can be obtained by using methods similar to those used in the classical 
case (Heyvaerts & Norman 2002a), giving: 



pr 



d a da d a da 
dz pr dz dr pr dr 



n 2 



p pc z 



d 2 a + 1 d ^da + fi '| Va | 2 
dz 2 r dr dr Q 



lE\a) 



Q'{a)p^ 



d Ida d Ida 

dz r dz dr r dr 



T up Va • V(<2p r_1 ) 



+ 



a' poa 2 p ( L — 7^r 2 f2 



a 



P 



[L' - 7 £rV) 



r - 1 r-ic 2 

L-^r 2 tt\ LL' 



|Va| 



(23) 



p a 2 £-pj r 2 £ " \p Q a 2 £-p) r 2 £ 

We find that, in comparison of the non- relativistic transfield equation, Eq.(23) has an ad- 
ditional second order term proportional to (Q 2 /c 2 ) on its left hand side which represents 
the cross-field component of the electric force. The gas pressure terms on the right hand 
side differ slightly from those obtained by Li (1993), who defined the function £ as being 
(1 + f dP/pc 2 ). We prefered the definition (5) which coincides with that of Li when the 
entropy Q(a) is independent of a. Toroidal variables, still implicit in 7, can be eliminated 
entirely by using the expression for the Lorentz factor in terms of poloidal variables (see 
section 2.3). The curvature of poloidal field lines is given by: 



Using the relation 



it- V)f = n d 4- 
as 



(V x Bp) ■ e e = n- V\B P \ - \B P \ d^/ds. 



(24) 



(25) 



- 9- 



the following equation, equivalent to Eq.(23), is obtained: 

K |Va| KK R } ? 7 r |Va| <9r pr 2 |Va| v 2/i c 2 2p 7 1 ^ 

The forces associated with the terms on the right hand side are the same as for classical dy- 
namics, except for the very last one which is part of the projection of the electric force normal 
to the magnetic surface. Another part of this electric force appears as a term proportional 
to Vt 2 r 2 /c 2 on the left hand side of Eq.(26). 

2.5. Force-free Limit 

The force-free limit applies in the case of very large magnetization. It corresponds to 
a limit in which the inertia p approaches zero. By Eq.(7) this implies that a approaches 
zero such that a/p remains finite. In this limit £, given by Eq.(5) reduces to unity and 
(po« 2 £/p) <C 1. Eqs.(19),(9),(10) then imply that the energy is all in Poynting form. Because 
the ordering between p and the critical density pqq 2 is opposite in the asymptotic limit, 
the latter is entirely out of the scope of the force-free approximation. The asymptotic limit 
applies to a state of the flow reached in regions much further away from the wind source than 
the limit down to which the force free approximation is valid. Nevertheless, the asymptotic 
regime is such that the component of the Lorentz force perpendicular to magnetic surfaces 
vanishes, except at boundary layers. There is no contradiction here (Nitta 1995). This 
asymptotic property occurs because the least negligible asymptotic forces normal to the field 
are electromagnetic. This property does not result from any a priori assumed dominance 
of electromagnetic forces over all the other forces present. The vanishing of the normal 
component of the electromagnetic force simply results from the wind dynamics. Between any 
near-source region in which force free conditions apply (because of strong magnetization) and 
the asymptotic domain (in which the component of the electromagnetic force perpendicular 
to the field vanishes as a consequence of the dynamics), an intermediate non- force- free region 
must exist. In this intermediate non- force- free region, currents cross magnetic surfaces and 
plasma is accelerated. The conservation of poloidal current on a magnetic surface (which, 
under force-free conditions, apply near the wind source) do not retain validity continuously 
out to the asymptotic domain. 

It is instructive to analyze how the inertia-less limit turns the transfield equation (23) 
into the well known force-free relativistic wind equation. In the inertia-less limit, Eq.(9), 
multiplied by a, reduces to — rB e = p aL = p I(a), where, as a approaches zero, aL 
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approaches the finite limit 1(a) which means that the poloidal current follows magnetic 
surfaces in the force-free regions and the associated Lorentz force has no component along 
the magnetic field. Multiplying Eq.(10) by a and taking the limit of vanishing a leads 
to —tVLBq = p aE = P, where it is meant that p^aE approaches a finite limit P as a 
approaches zero. This implies that the energy flux is all in Poynting form. This relation can 
be restated as aE = IQ. The point r defined by Eq.(18) reduces in this case to the light 
cylinder radius r — c/Q. In the force-free limit, Eq.(23) should be expanded in p and a to an 
appropriate order, to take account of the cancellation of the dominant terms. The inertial 
terms on the left hand side of Eq.(23) become negligible. The pressure terms disappear from 
its right hand side and aL approximately equals I. One is eventually left with the well 
known force-free pulsar wind equation (Beskin et al. 1993): 

/ OV\ (d 2 a d 2 a\ ( n 2 r 2 \ 1 da rM' l2 2 TTf/ s 
( 1 " — + ^H 1 + ^ r ^ ~ ^ |V01 + ^ 11 (0) = ° (27) 



3. Field Regions 

3.1. Transfield Equation 

Let us compare the terms on the right of Eq.(23) in the large-r limit. The gravity term 
and the centrifugal force term, which declines as 1/r 3 , become negligible, and the pressure 
term becomes very small, so that £ approaches unity. The electric force remains of the same 
order as the hoop stress, though. The asymptotic form of the relativistic transfield equation 
is then: 

2 2 / pr 2 n 2 \ # lVo f~f Vo 2 ^ r \ fr 2 B 2 Q 2 r 2 \Va\ 2 \\ ^ 

^ { 1+ ^)i-pW\V[^ + Qp ) + ^ " s-^) ) (28) 

Eq.(10) shows that, on a given magnetic surface, tBq is bounded at large distances and 
Eqs.(12) and (7) show that pr 2 and r\Va\ are also bounded. When they approach finite 
limits, the hoop stress term in Eq.(28) decreases as 1/r, as does the electric force. Both the 
poloidal magnetic pressure and the gas pressure decrease more rapidly with r. It has been 
shown (Heyvaerts & Norman 1989) that the inertia force associated with the curvature of 
the poloidal motion on the left of Eq.(28) must decrease faster than 1/r. It becomes negli- 
gible compared to the hoop stress and electric force. The asymptotic form of the transfield 
equation becomes: 



Va-V[r 2 B 2 -r 2 B 2 P ^) =0 (29) 



Eq.(29) generalizes to the relativistic case our earlier non-relativistic result (Heyvaerts & 
Norman 1989, 2002a; Chiueh et al. 1991). Adding to the left hand side of Eq.(29) that 
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negligible part of the electric force which is proportional to the curvature of poloidal field 
lines the relation, we find: 

Va- (jx B + p e E^j = (30) 

So, the component of the electromagnetic force normal to the magnetic surface asymptot- 
ically vanishes on flared surfaces. This does not imply a strictly force-free situation since 
Eq. (30) is asymptotically satisfied by the vanishing of each of its terms separately and holds 
only perpendicular to field lines. We refer to regions where Eq.(29) holds true as field-regions 
of the asymptotic domain. The pressure force becomes significant near the polar axis and 
near neutral magnetic surfaces, where the electromagnetic force vanishes. Eq.(29) can be 
integrated following orthogonal trajectories to magnetic surfaces. Let b label one such or- 
thogonal trajectory. We define b as being the distance to the origin of the point where this 
orthogonal trajectory meets the polar axis. The integrated form of Eq.(29) is: 

r 2 B 2 9 -^\Va\ 2 = plK 1 (b) (31) 

K\ is independent of a. Using Eq.(ll) and noting that Be ^> Bp in the large-r limit, Eq.(31) 
can also be written as: 

r 2 (c 2 B 2 -E 2 ) =^c 2 K^b) (32) 

The scalar invariant of the electromagnetic field (c 2 B 2 — E 2 ) appears on the left of this 
equation. From the asymptotic form of Eqs.(7) and (12), the toroidal field can be expressed 
in terms of poloidal variables as: 

7oopr 2 fi a Qr\Va\ 

rB e = = - — (33) 

a \a\ 

When used in Eq.(31) this gives: 

p 2 K 1 (b)=r 2 B 2 (l-^) (34) 

Then, Ki(b) is positive and can be written as: 

K^b) = K 2 {b) (35) 

K{b) has the dimension of an electric current. In the non-relativistic case the quantity which 
becomes independent of a in field-regions is the total poloidal electric current I. Eq.(34) 
indicates that the situation is different in a relativistic plasma flow. K and / are related by 
(Chiueh et al. (1991)): 

I{a,b)= loo {a,b)K{b) (36) 
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This defines K as an algebraic quantity having the sign of I (see Eqs.(14) and (33)). Since 
the azimuthal velocity asymptotically approaches zero, the terminal Lorentz factor, 700(a), 
refers to the terminal poloidal velocity v^a). The Bernoulli equation (10) becomes in the 
same limit: 

7 ooc 2 = E- — (37) 
a 

Eqs.(36) and (37) provide expressions for the current / and the momentum 700^00 in terms 
of K(b) and the first integrals: 

'M^' ^SiH (38) 



700 00 (c' + KOIa) { ' 

K is the poloidal current observed in a rest frame where the poloidal motion vanishes. This 
can be seen by transforming from the laboratory frame to a frame moving with the fluid at 
the poloidal velocity Voo along the direction of the poloidal field. The azimuthal magnetic 
field is given by Eq.(33) and the electric field by Eq.(ll). K is equal to: 

/io7oo I a I ^o7oo^oo 

In the moving fluid rest frame, the electric field vanishes, while the azimuthal magnetic field 
is: 

PoK 

Bejiuid — -B6>/7oo = — (41) 

The arc length element (r dO) remains invariant in the transformation, so that: 

r d6 Bejimd = -p Kd9 (42) 

This shows that K is negatively proportional to the current enclosed by a circle of radius r 
carried by the fluid motion. We refer to K as the proper current. 



3.2. Current-Carrying Boundary Layers and Electric Circuit 

— * 

Eq.(28) shows that the electromagnetic force is proportional to KVK. Since K vanishes 
with /, this force yields to pressure at boundary layers around null surfaces and near the 
polar axis. However, since the pressure is weak, the thickness of the boundary layers must 
be small. Using Eq.(33), Eq.(28) reduces to: 

I(V«.V) (Qff) = -±(Va- V) C'^l - pi IVaf) = - (^) (43) 

p pr 2 \ 2p 2p c 2 J a \A«o7oo*W 
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The proper current K then has the following profile along an orthogonal trajectory: from 
zero at the polar axis, it quickly rises to a non-zero value at the edge of the circum-polar 
boundary layer, then stays constant and steeply returns to zero through a boundary layer 
about the next null surface. The current system closes exactly in cells bordered by null 
surfaces. Eq.(106) of section (6.3) shows that K 2 resumes its original value after crossing 
the boundary layer about a null surface. K only changes sign. In the field-region of the next 
cell, K remains constant, again returning quickly to zero at the next null surface. 



3.3. Asymptotic Grad Shafranov Equation 

Eq.(29) can be transformed to give: 

d ( a da\ d ( a da\ _ -* ( alVal \ alVal dip 



dz \pr dz J dr \pr dr J \ pr J pr ds 

Since in the asymptotic domain (r dip /ds) becomes negligibly small, the term on the right of 
Eq.(44) can be neglected compared to any one of those on its left. Expanding second order 
operators and using Eq.(7), this leads, similarly to the non-relativistic case, to the following 
equation for a(r,z): 

Aa = Va- V(ln(r|Va|)) (45) 

which can also be restated as: 

div 0) = (46) 
Multiplying Eq.(45) by r 2 and denoting: 

Va - V d 



|Va| 2 da 



(47) 



it can also be transformed into: 



da V 2 

Using Eqs.(31), (33) and (35), this can finally be restated as a Grad Shafranov equation: 

The boundary conditions to Eq.(49) are that a = along the polar axis and that a = A on 
the equatorial plane. When K(b) is constant and non-zero, these conditions imply that a 
depends on the latitudinal angle ip only. In this case, Eq.(49) becomes an ordinary differential 
equation for a(ip) which reduces to the form of Eq.(55) below. 
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4. Solution in Field Regions 
4.1. Hamilton Jacobi Equation 

Using Eqs.(33) and (39), Eq.(31) can be restated as: 



r Va = rv \( 2 i YFTTmTi \ i ( w = f{ a , b ) (50) 

Let a be the curvilinear abscissa along an orthogonal trajectory to magnetic surfaces, con- 
ventionnally increasing from pole to equator. An equivalent form of Eq.(50) is: 

da _ (c 2 + K(b)n(a)/a(a)) tt(a)da 

~ ~ fi \K(b)\y/E^a) - (c 2 + K(b)Q(a)/a^f 

When K(b) approaches a finite constant at large distances, Eq.(50) can be further 
transformed, by defining S(a) = J^da'/f(a'), into the following Hamilton- Jacobi equation: 

\VS\ = - (52) 
r 

the solution of which can be constructed by ray tracing, as in the non-relativistic case 
(Heyvaerts & Norman 2002a). The boundary conditions at the equator have been shown to 
select a solution in which orthogonal trajectories to magnetic surfaces are circles centered at 
the origin. This is only an asymptotic, approximate, result, as is Eq.(52) itself. 



4.2. WKBJ Approximation 

When K(b) declines to zero at large distances, the wind becomes asymptotically kinetic- 
energy-dominated. The function f(a,b) of Eq.(50) asymptotically vanishes in this case. 
Eq.(50) then does not give r|Va| as a function of a only. If, however, the decline of K(b) with 
distance is very slow, the WKBJ approximation allows to neglect this variation. Orthogonal 
trajectories then remain quasi-circular. In the vicinity of the orthogonal trajectory of radius 
b, the flux surfaces are approximated by a series of nested conical surfaces locally represented 
by 

z = r tan(tp(a,b)) (53) 
The angle ip(a, b) is assumed to slowly vary with b. 
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4.3. Solution in Field Regions 

Eq.(50) is now considered in the WKBJ approximation and in the geometry described 
by Eq.(53). We find: 

m i cosib 
r \aip/da\ 

This gives the following differential equation for ip at given b: 

dip n(a) da (c 2 + K(b)fl(a) /a(a)) 



costp fioK{b) c ^E 2 {a) - (c 2 + K{b)Q{a) /a{a))< 
which integrates to: 



(55) 



/ x / ,n , ( f ai 1 n(a') da' (c 2 + KQ(a')/a(a')) \ . . 

tan^ (a, b) = tan^ ai, 6) + sinh / — / ; v y !l v =U= 56 

\Ja c y/E 2 (a>) - (c 2 + KQ(a')/a(a')y J 

where a\ is a reference flux in the cell under consideration and, again, K depends weakly on 
b. If the cell begins at the equator, ai can be taken as the flux variable, A, of the equatorial 
surface and t&n(ip(A, b)) = 0. This neglects the small flux in the equatorial boundary layer, 
if the latter is a null surface. These results are similar to those obtained for non-relativistic 
winds by Heyvaerts & Norman (2002a) and for relativistic winds by Nitta (1995). 



4.4. Flux Distribution in Cylindrical Regions of the Field 

When K(b) approaches a finite limit K^, there may exist regions of the free-field where 
magnetic surfaces become cylindrical. Their radius is given by Eq.(51), which in this geom- 
etry gives: 



/ r n(aVaV + A^(a-)/a M ) | 



where a2 is a reference flux in the cylindrical field-region. Eqs.(56) and (57) give approxi- 
mately identical results when tan^ becomes very large. 

At this point the solution in the free field is described by Eq.(56). To make this solution 
complete, we must solve Eq.(43) in boundary layers and determine how the circumpolar 
proper current K(b) depends on b. 
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5. The Polar Boundary Layer 

5.1. Solution in the Polar Boundary Layer 

Plasma pressure, or possibly poloidal magnetic pressure, must be taken into account 
in the vicinity of the polar axis. Since the poloidal magnetic pressure decreases faster with 
increasing r than the plasma pressure, the latter dominates, unless the polytropic entropy Q 
vanishes. The transfield force balance near the pole is between the hoop stress, the electric 
force and the pressure gradient. Using Eq.(33), Eq.(43) takes the form: 

(Va • V) (^QP T - 1 ) + -(Va • V) (e™) = (58) 
VI - 1 J a V Po a J 

Assume that Q, Q and a can be taken as constants in the boundary layer. This is valid for 
small K (Heyvaerts & Norman 2002a). Eq.(58) then integrates to: 

' -Qo/- 1 + ^7 = K M) (59) 



r - 1 /io«o 

The integration constant K 2 (b) can be identified by considering the left hand side of Eq.(59) 
on the polar axis, where the density is po(b), so that 

K 2 (b) = J—Q oP l-\b) (60) 

Alternatively K 2 (b) can be identified by considering the left of Eq.(59) at large axial dis- 
tances, where the pressure becomes negligible. Using Eqs.(l), (7), (33) and (31) this gives: 

K 2 (b) = ^K(b) (61) 

where K(b) is the proper current, defined by (Eq.(35). Eq.(59) can be solved for r in terms 
of the parameter 

x = p/p (b) (62) 

This results in Eq.(63). Using this expression for ip in terms of x in Eq.(54) to express a as 
a function of x, a parametric representation of the flux distribution in the polar boundary 
layer is obtained: 
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The Lorentz factor and velocity 70(6) and vo(b) refer to their values at the polar axis at a 
distance b from the source. Eqs.(lO) and (5) give 70 as: 

^ [ c2 + f^Qop r o~ 1 (b)^=E (65) 

from which we get 

lovo = c (66) 

c 2 + — Q0P0 (b) 



5.2. Matching the Polar Boundary Layer Solution to the Outer Solution 

The inner limit (a <C A) of the outer solution (Eq.(56)) must now be asymptotically 
matched to the outer limit (x <C 1) of the inner solution, Eqs.(63)-(64). These equations 
provide , for x <C 1, the following expression for if) (a, b): 

2,, H r Qo gottg p r - 2 ( 2a (T -1)0.1 (f+jzjQvPt*) \ 

» s 6 » - ( r - 1) W «» (ro^^ c 71 - (g + g QoP ;-V j ' ' 

where po depends on 6. In the vicinity of the polar axis, tan^ is large. For very small cosijj, 
the first term on the right of (Eq. (56)) is negligible and this relation can be written as: 

costp = — - — (68) 
coshx 

where the hyperbolic argument x( a > b) is: 

y[a b ) = 1 f 1 (c 2 + K(b){l(a')/a(a>))n(a>) da' 
fi cK{b) J a y/E 2 (a') - (c 2 + K(b)Q(a')/a(a')) 2 

The apparent weak dependence of x 011 a i is absorbed by the neglected first term on the 
right of Eq.(56). The relation (71) is thus essentially independent of a±. For definiteness, 
ai can be taken, in the case of dipolar symmetry, as the equatorial value, A, of a. For very 
small a: 

(ab , 1 ( f ai (c 2 + Kn(a')/a(a')) Qjaf) da' (c 2 + ROp/ap) n a , 

X[a ' ~ fi cK yj ^E 2 (a')-{c 2 + KQ{a')/a{a')) 2 ^E 2 - (c 2 + KQ /a ) 2 1 

From Eqs.(68) and (70) the inner limit of the outer solution can be written: 

cos 2 V; =ex / 2_f f ai (c 2 + KQ/a) Qda' _ a Q (c 2 + KQ /a ) > 

4 6XP ^ I2 cK\J v /E 2 -(c 2 + KQ/a) 2 " ^ E 2 - (c 2 + KQ /a ) 2 ] 1 

For brevity the dependence of K on b and of the first integrals on the integration variable a! 
has been omitted. 
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5.3. Bennet Pinch Relation 

For the solutions (67) and (71) to smoothly match, it is necessary that the arguments 
of their exponential functions of a coincide. This implies: 

-L-Q fi-\b) = K(b) ^ (72) 
1—1 cto 

This is again Eqs.(60)-(61). Eq.(72) is a Bennet pinch relation between the proper current 
and the plasma pressure on the axis. The proper current K is related by Eq.(36) to the 
asymptotic poloidal current and Lorentz factor. 



5.4. Polar Boundary Layer Proper Current, Density and Radius 

For a smooth asymptotic matching of Eqs.(67) and (71) the factors in front of their 
exponential functions of a must also coincide. Let us note 

s " = (f3i) G°Pb -1 (73) 

The velocity s is of order of the sound speed at the axis. It depends on po(b). Taking Eq.(72) 
into account, this matching condition can be written as: 

1 Q Wo = 4 exp / U (74) 

(r - 1) Vt 2 b 2 \ c/i a s 2 Jo A !E 2 (a') (c 2 I d 1 "Kfo ^ I 



a(a')Qo - 



Eq.(74) determines p (b) (on which s 2 depends by Eq.(73)). Let us define a length £, a 
dimensionless measure of the axial density, n , and a reference magnetic flux a by: 

P 2 _ T g (/ioag) r ~ 1 

1 ~ r=i % (75j 

n {b) = g (76) 

ao = ^ c poa f (77) 



Let us also introduce the notation 



„2 , c .2 n(a')a 



Jo c ^E 2 {a>) - {c 2 + s 2 ggg) 2 fi oa 
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The integral A depends on n (b), since, by Eq.(73), s 2 does. The logarithm of Eq.(74) 
provides the following equation for n (b): 

M«o(&))_, ^„ ^ , , f& 



(2-r)ln(no(6))+ln^ (79) 



Since n (6) is small and (b/£) large, a solution can be obtained by iteration, giving, in the 
simplest degree of approximation: 



The solutions of Eqs.(79) or (80) are controlled by the growth of the logarithmic term on 
their right. They can be satisfied for large 6's in two different ways, according to whether 
the proper current K(b) approaches a finite value or decreases to zero. If K(b) approaches 
a finite value, Eq.(72) indicates that the axial density becomes independent of distance. 
The logarithmic term in the denominator of equation (80) must then be compensated by a 
divergence of the numerator term, X(n (b)). As b increases, K(b), related to n (b) by Eq.(72, 
should approach a limit that causes the integral A(n ) to diverge. This implies that K should 
approach the flat absolute minimum, K sup , of the function (a(E — c 2 )/Q). If this function 
does not have such a minimum, a solution with a finite limit for K(b) cannot be found and 
K(b) should actually approach zero. Because the square root at the right of Eq.(78) must 
remain definite for any a, K sup can only be approached from below. When K = K sup , all 
the energy flux is in Poynting form. 

If K(b) declines to zero at large distances, A approaches a limit Ao independent of b. 

i 

Eq.(80) then shows that p (b) scales as (\n(b/£)) r - 1 and K(b), given by Eq.(72), slowly 
decreases as (hi(b / : 

m = ( l A M i (si) 

This slow decrease of K(b) justifies a-posteriori the adopted WKBJ procedure. The hyper- 
bolic argument x °f Eq.(69) becomes, in this case (taking a\ equal to the equatorial value of 
the flux, A): 

where 



Ja 



fc(fl) = W*pEt (83) 

L A da> n ^v^ 
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The solution (63) shows that the core radius of the jet at a distance b from the source is: 

rl(b) = - JL Q^py^ (84) 

At this point the solution near the pole and in the field region extending in the polar-most 
cell is completely determined. 

It is important to note the very slow decline of the proper current and associated Poynt- 
ing flux that Eq.(81) implies. The ratio a of the Poynting flux to the other forms of energy 
flux is given by Eq.(10): 

a = rQ(a)B e 

/loa (E(a) + r£l(a)Bg/ pioa) 

From Eqs.(14) and (38) it results that 

ff= ^w (86) 

Eq.(81) then explicitly indicates how a decreases in a kinetic wind with distance b to the 
wind source, i.e.: 



/i c a(a) [Jo ^n 0V ^V)^j W(b/i) 




(87) 



the convergence of the flow to a completely kinetic energy dominated state is obviously 
very slow. Chiueh et al. (1998) stress the fact that such a slow decline is a difficulty for 
understanding the Crab pulsar wind, which has a very large terminal Lorentz factor and a 
cr-parameter as low as 10~ 2 - 10 -3 . We come back in paper III of this series (Heyvaerts & 
Norman 2002b) on the implications of such a slow decline and show that, even though the 
mathematical asymptotics implies a very low a, actual winds may not have converged to this 
state during their finite life time, so that a could remain finite when they reach, for example, 
a terminating shock. We do not propose in this paper any specific solution to the question 
raised by Chiueh et al. (1998) concerning the properties of the Crab pulsar wind. Such an 
explanation may have to be found, as Chiueh et al. (1998) suggest, in non-ideal effects close 
to the light cylinder not considered here. 



5.5. Force- Free Polar Boundary Layers 



In the exceptional case when the poloidal magnetic pressure dominates over plasma 
pressure in the asymptotic circumpolar region, the transfield mechanical balance equation 
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takes the form of Eq.(28), neglecting the poloidal curvature term and pressure. In cylindrical 
geometry, this gives (Appl and Camenzind 1993a,b): 

dB 2 1 d ( 2 „o fi 2 r 4 2 \ . . 

Under negligible gas pressure, £ = 1. Eqs.(7), (10) and (12) provide B P in terms B e : 

c 2 B 2 / c 4 



Bp r 2 fl 2 i, 1 (E + rQBe/pioa) 2 ) (§9) 
Eq.(88) can then be restated as an equation for I defined by Eq.(14): 
d ( I 2 ( c 4 \\ \ d ( I 2 c 2 



dr Vr 4 ^ 2 V - /^/a) 2 J J r 2 dr \(E - IQ/a) 

Assuming the first integrals to be constant in the polar boundary layer, which implies that 
I <^aE/fl (Heyvaerts & Norman 2002a), Eq.(90) simplifies to: 



E 2 -c A d fl 2 \ 1 dl 2 



c 2 VLq dr \ r 4 / r 2 dr 



the solution of which is: 

^^ w+IS-c-) (92) 



The characteristic thickness r of the axial pinch is then: 

2 % 



r °-w- (93) 



Eq.(92) establishes a relation between the total electric current 1 00(b) flowing through the 
polar boundary layer and the limit of (I/r 2 ) as follows: 

For small 1 00(b) the density and magnetic field on the axis are given by: 



When the first integrals are almost constant in the boundary layer, K^b) is, from Eq.36: 

*«,(&) = %^ (97) 
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Using Eqs.(95), (96) and (97), Eq.(94) reduces to: 

BUb) = noK^b) (9g) 
li p (b) a 

For jets having a force-free polar boundary layer, this equation replaces Eq.(72). Compared 
to Eq.(72), which applies to gas-pressure-supported polar boundary layers, the axial Alfven 
velocity has replaced in Eq.(98) the axial sound speed. 



6. Null-Surface Boundary Layers 

6.1. Divergence of Mass-to-Magnetic Flux Ratio at Null Surfaces 

At a null surface of flux parameter a n , the mass flux to magnetic flux ratio, a (a), defined 
by Eq.(7), diverges if there is a mass flux on this null surface. It has been shown in Heyvaerts 
& Norman (2002a) that near a n the function a (a) behaves as 

a(a) « -. ] —— (99) 

w \a n - a\ v v ' 

where v is positive and stricly less than unity. In the case, assumed below, of a field of a 
dipolar type of symmetry, the null surfaces reduce to only the equatorial plane. We now 
outline the solution near the equator. 



6.2. Solution in the Equatorial Boundary Layer 

The structure of the flow is given by Eq.(58). We assume that the first integrals Q and 
Q are almost constant in the boundary layer, with values Q e and Q e but we cannot assume 
this for a, because of its divergence. Taking into account the vanishing of (Va • Vr) at the 
equator, Eq.(58) becomes: 

*-*K + S) = (100) 

R enters Eq.(lOO) as a parameter since (Va • VR) vanishes. Eq.(lOO) integrates to: 

Q e p r + = QePe(R) (ioi) 

p e {R) being the equatorial density at the distance R. The distribution of magnetic flux in 
the equatorial sheet can be found in terms of the parameter 

x ^ m (102) 
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From Eq.(lOl), a can be expressed at a given R as a function of X by: 

1 2Q e p r e - 2 (R) ( 1 1 



p a 2 (a) tt 2 e R 2 \X 2 X 2 ~ r 



(103) 



This implicitly determines a(X) since a is supposedly known as a function a. The distribu- 
tion of flux with latitude angle ip, at almost constant R, is infered using Eq.(7). Eqs.(54), 
(103) and (10) yield the following differential equation for ip at given R: 



Although not given here, ip(X) can be determined by quadratures by using Eq.(103) which, 
for known, and locally invertible, a(a) provides a as a function of X. 



6.3. Matching the Equatorial Boundary Layer Solution to the Field 

In the equatorial boundary layer, the solution is expressed by Eqs.(102), (103) and (104). 
This boundary layer solution must asymptotically match the field-region solution, expressed 
in differential form by equation (55). The outer regions of the boundary layer correspond to 
small X. In this limit, X can be eliminated by using Eq.(103), so that Eq.(104) reduces to: 

= „ II. da ^+j^SSS (105) 

Matching requires that equations (55) and (105) become identical, implying 

fx K(R) = ^-J2^Q e R 2 pl(R) (106) 
\a\ 

This expresses the balance between the gas pressure force and the electromagnetic force 
across the equatorial boundary layer. When the proper current, K(R), approaches a finite 
limit as R grows to infinity, the equatorial density decreases as 

Pe (R) ~ R-t (107) 

For kinetic winds, K(R) decreases as (ln(i?/£)) _1 . The equatorial density then declines with 
distance as 

p e (R) ~ (R HR/t)Y$ (108) 
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7. Shape of the Magnetic Surfaces 



We have now obtained a complete solution in the asymptotic domain, both in field 
regions, near the pole and near null magnetic surfaces. The dependence of K(b) on b is 
obtained from Eqs.(72) and (79). We now possess all the information needed to calculate 
the shape of magnetic surfaces in any region of the asymptotic domain. 



The magnetic surfaces of Poynting jets near the polar axis are cylinders. Their radius 
is given by Eq.(57) (with Eqs.(84) and (72)) if they are in the free field, and by equations 
(63)-(64) if they are in the polar boundary layer. The outer-most magnetic surfaces in such 
flows are conical. The angle ip defined by Eq.(53) becomes independent of b and is given 
by Eq.(56). How cylindrical and conical regions smoothly merge one into the other is to 
be discussed in a forthcoming paper (Heyvaerts & Norman 2002b) . The shape of magnetic 
surfaces in the equatorial boundary layer of Poynting jets is discussed in section (7.4). 

7.2. Magnetic Surfaces in Field Regions of Kinetic Winds 

The shape of poloidal field lines in kinetic winds is described by the differential system 



with ip(a,b) given by Eq.(56), supplemented by Eqs. (81) or (82). In the case where tp is 
close to it/ 2 we obtain: 



where k(a) is defined by Eq.(83). The magnetic surfaces are a collection of nested power-law 
paraboloids of variable exponent q(a) = (1 — k(a)) as represented in Figure 1. The shape 
of magnetic surfaces for which the approximation ip m it/ 2 is invalid can be treated as in 
Heyvaerts & Norman (2002a), with a similar result. 



7.1. Magnetic Surfaces in Field Regions of Poynting Jets 



dr = cosip(a, b) db 
dz = sin if) (a, b) db 



(109) 




(110) 
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7.3. Magnetic Surfaces in the Polar Boundary Layer 

The solution deep inside the polar boundary layer of kinetic winds is obtained from 
Eqs.(63)-(64) in the limit of x close to unity. The axial density, po(b) is given by Eq.(79) 
with the factor A now being equal to 

A = ! A n{a ' )c2 ^ fill) 
° Jo n ^E^(a')-c 4 a 



Eq.(63) can be written as 

r 2 1 



i 2 n 2 - r (b) \x x 2 ^ ^ 
while Eq.(64) provides the value of the parameter x in terms of a by: 



(112) 



,l x 2-T, r ,A a c 2 1 



r-1 J ao V E o - c 4 nl'^b) 

From Eq.(80), the dimensionless axial density is in this case given by: 

1 In (4b 2 /I 2 ) 



nl-\b) A 



(114) 



When x is close to unity, it can be eliminated between Eqs.(112) and (113), taking b to be 
almost equal to z. This gives the shape of magnetic surfaces in this region as: 



Z V°o (£ 2 -c 4 )3 V A o 



' " f ^ 2V(r " 1} f m 2 4)Y T - 1) (lis) 



i 



7.4. Magnetic Surfaces in the Equatorial Boundary Layer 

The information on the shape of magnetic surfaces in the equatorial boundary layer 
is contained in the parametric solution provided by Eqs.(103)-(104), the equatorial density 
p e (R) being related to K(R) by Eq.(106). The latter approaches a non-vanishing constant 
for Poynting jets and decreases logarithmically for kinetic winds. It can therefore be written 

as 

K ^ - www (116) 

where m = for Poynting jets and m — 1 for kinetic winds. The factor J m is different in the 
two cases. Eq.(106) then gives the associated equatorial density. For small X, Eqs.(103)- 
(104) give the results for the field-region. In this limit, a can be explicitly obtained from 
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Eq.(103) in terms of X and Eq.(105) integrates to: 

(117) 



z 



A a Q e da S + j^^QJPr/JR) 



n J. H VOTVFW c ^J^==y 

For Poynting jets, K(R) approaches a finite constant as does R 2 p T e (R) (from Eq. (106)). Con- 
sequently, the magnetic surfaces become conical at the outskirts of the equatorial boundary 
layer. For kinetic winds K(R) decreases logarithmically and the equatorial density declines 
as (R ln(R/£))~T . Eq.(117) shows that in the equatorial boundary layer z is proportional to 
(A — a) R\n(R/£). The magnetic surfaces at the outskirts of the equatorial boundary layer 
become slightly convex paraboloids. By contrast, deep inside the equatorial boundary layer, 
gas pressure dominates, X is close to unity and Eqs.(103) and (104) give: 

, = r ^ , (118) 



In this region z(R) scales, at fixed a, as z ~ R~^~ for Poynting jets and as z ~ R~^~ (in-j) r 
for kinetic winds. Magnetic surfaces are concave, i.e. they bend towards the equator. This 
result is quite consistent with the properties of such flows discussed by Okamoto (1999, 2000, 
2001) and Beskin & Okamoto (2000). It does not contradict our earlier results that the 
poloidal field lines cannot asymptotically bend towards the equator (Heyvaerts & Norman 
1989). Actually, any such poloidal line eventually escapes the equatorial boundary layer 
region at a finite distance, rejoining the field region where it becomes convex, bending away 
from the equator. It can be shown that when poloidal field lines exit the equatorial boundary 
layer, they do so at an angle to the equatorial plane that decreases with distance. The proof 
is similar to the non-relativistic case (Heyvaerts & Norman 2002a). 



8. Conclusions 

Our main results are summarized as follows: 

I. We have derived a global solution for the asymptotic structure of relativistic, stationary, 
axisymmetric, polytropic, unconfined, perfect MHD winds, assuming their five lagrangian 
first integrals to be known as a function of the flux parameter, a. Current-carrying boundary 
layers along the polar axis and at null magnetic surfaces are features of this solution, which is 
given in the form of matched asymptotic solutions separately valid inside and outside these 
boundary layers. 
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The asymptotic structure consists of field regions where, although the electromagnetic force 
decreases to zero asymptotically, it still dominates the force balance. In these field-regions 
the magnetic field structure is shown to be described by the Hamilton- Jacobi equation (52) 
which we solved. We obtained the distribution of flux as represented by equations (56) and 
(57). The poloidal proper current K (Eq.(40)) is shown to remain constant in each of these 
field-region cells, approaching a value independent of a. For a non- vanishing asymptotic value 
of K a substantial part of the energy reaches the asymptotic domain in electromagnetic form. 
For vanishing K, all of the wind energy is in kinetic form. We found that K may approach 
zero as the inverse logarithm of distance. 

Field regions are bordered by 

1. boundary layers along the polar axis and 

2. null surface boundary layers 

These boundary layers have the structure of charged column or sheet pinches supported by 
plasma pressure. For very cold winds, the support is by magnetic poloidal pressure. They 
carry the poloidal electric current still present in the asymptotic domain. 

II. The polar boundary layer has been shown to have the structure of a charged column- 
pinch for which a Bennet relation between axial pressure and proper axial current is given by 
Eq.(72). The distribution of flux with radius in this region is explicitly given by Eqs.(63)- 
(64). 

III. The equatorial region, or more generally the vicinity of any null magnetic surface, has 
been shown in section (6) to have the structure of a charged sheet-pinch. The density in the 
equatorial boundary layer of a Poynting jet decreases as described by Eq.(107), while in the 
equatorial boundary layer of a kinetic wind it decreases as described by Eq.(108). 

IV. We have derived solutions valid in these regions and have matched them asymptotically 
to the field region solutions. Matching between the polar boundary layer and the field-regions 
determines the decline with distance of the density at the polar axis or, equivalently, the 
proper current around the polar axis. The variation with distance of the axial density is 
given by the solution of Eq.(79), considering Eqs.(74),(75), (76), (77) and (78). A simplified 
form of Eq.(79) is Eq.(80). 

V. The geometry of magnetic surfaces in all parts of the asymptotic domain has been explic- 
itly deduced in terms of the first-integrals in section (7) for both kinetic winds and Poynting 
jets regimes. In particular, the magnetic surfaces in field-regions of a relativistic kinetic wind 
are paraboloids with a variable exponent, given by Eqs.(llO) and (83), while, deep in the po- 
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lar boundary layer, magnetic surfaces have exponential shapes, as represented by Eq.(115). 
Deep in the equatorial boundary layer, the shape of magnetic surfaces is represented by 
equation (118). Note that any magnetic surface embedded in this equatorial boundary layer 
eventually finds its way out at some finite distance. 

These results are illustrated in Figures 1-7. Figures 1-3 represent properties of kinetic 
winds as deduced from our analysis. Figure 1 illustrates the geometry of magnetic surfaces. 
Figures 2 and 3 represent the structure of the polar and equatorial boundary layers. Figures 
4-7 represent similar properties for Poynting jets. 

In all cases, the polar and null surface boundary layers which carry residual electric 
current may stand out observationally against the field-regions, both because of their large 
density contrast to them and because they possess a source of free energy that makes them 
potentially active by the development of instabilities. It is noteworthy that the density 
about the pole either does not decline with distance in the case of Poynting jets, or only 
very slowly (as a negative power of the logarithm of the distance) in the case of kinetic 
winds. It may be that what is observed as a jet be but the dense and active polar boundary 
layer of a flow developed on a much larger angular scale. This diffuse flow may be itself 
barely visible because of its very low density and activity. Null surface boundary layers, for 
example equatorial ones, enjoy a more favourable status than field regions from this point 
of view (Beskin & Okamoto 2000), and may be observed in association with jets. The X-ray 
structure of the Crab Nebula (Weisskopf et al. 2000) can be understood in such a framework 
(Blandford 2002). In addition, our work is relevant to the large scale aspects of pulsar winds, 
jets from active galaxies and gamma ray bursts (Vlahakis & Konigl 2003a,b). It is interesting 
to note that, from our analysis, the highest degrees of collimation are associated with flows 
which carry significant Poynting flux. In our next paper (Heyvaerts & Norman 2002b), we 
show that such flows should actually occur, owing to the slow decline of the asymptotic 
proper current with distance. 
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Fig. 1. — The magnetic field structure for kinetic winds in the asymptotic field region, from 
Eq.(llO). The functions E(a) and Q(a) have been taken as constant which implies that 
q(a) = (a/A). On the pole g(0) = and at the equator q(A) = 1. The field lines in each 
quadrant correspond to a/ A = 0.2, 0.4, 0.6, 0.8, 0.9. An interpolation formula has been used 
to connect the asymptotic solution to a split-monopole field at the origin. The scale for r 
and z is arbitrary. 
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Fig. 2. — The behaviour of kinetic winds near the polar axis. The panels show: magnetic 
field structure (central frame); the normalized density (Eq.(76)) at the polar axis vs. the 
distance along the polar axis (upper right); the total proper current integrated about the 
polar axis vs. the distance along the polar axis (lower right); the ratio of the density to 
its polar value at the same z across the polar pinch (upper left) and the integrated proper 
current across the polar axis (lower left). The latter two quantities are plotted for some 
arbitrarily chosen z . In the central frame, a = a (y) 3 ^ (Eq.(115)). The slightly parabolic 
shape of the field lines is not clearly visible on the scale of the plot. The right panel curves 
are from Eqs.(72), (75) - (77) and (80) with Ao given by Eq.(78) for negligibly small K^,. 
The left panel curves are from Eqs.(62)-(64), (33), (34), (75)-(77) and (80). 
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Fig. 3. — The behaviour of kinetic winds about the equatorial boundary layer. The panels 
show: the magnetic field structure (central); the normalized density n e = in the equa- 
torial plane vs. distance along the equatorial plane (upper right); the total proper current 
integrated about the equatorial plane vs. distance along the equatorial plane (lower right); 
the ratio of the density to its equatorial value at the same r across the equatorial sheet 
(upper right) and the integrated proper current across the equatorial sheet (lower left). The 
latter two quantities are plotted for some arbitrary value of r, r . For simplicity, Q e , a e 
and Q e have been set equal to Q , a and Q respectively. The field lines (central frame, 
Eq.(118)) correspond to ^ — 0.9. Their parabolic shape is not clearly visible on the scale of 
the plot. The plots in the right panel are from Eqs.(106) and (108), while those in the left 
panels are from Eqs.(99), (102), (103) and (33) with v = \. 
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Fig. 4. — The magnetic field structure for Poynting jets in the field region, from Eqs.(56) 
and (57). From pole to equator, the contours are for (a/A) = 0.4, 0.6, 0.7, 0.85, 0.9 and 
0.95. E(a) and Q(a) are constants, while we assume, as discussed in Paper III, that a(a) = 
c 2 )){l+ \((a/A)-^f/^l 



}. The unit of the plot for the variables r 
ation formula has been used to connect the 



(KQ/(E - c 2 )){l + [((a/A) - 0.8)7^1 - (a/A) 
and z is the scale £ defined in Eq.(75). An interpo 
asymptotic structure of the field to a split monopole field near the origin. Eqs.(56) and (57) 
become increasingly accurate with larger r/£ and z/£, the scale of the transition between the 
split monopole structure near the origin and the asymptotic one being arbitrary. 
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Fig. 6. — The behaviour of Poynting jets near the polar axis. The plots show: magnetic field 
structure (central frame); the normalized density at the polar axis vs. the distance along 
the polar axis (upper right); the total proper current integrated about the polar axis vs. the 
distance along the polar axis (lower right); the ratio of the density to its polar value across 
the polar pinch (upper left) and the integrated proper current across the polar axis (lower 
left). The latter two quantities are plotted for some arbitrarily chosen z . The central frame 
shows that the field lines are exact cylinders (Eq.(84)). The right panel curves use Eqs.(72) 
and (76). The left panel curves are obtained from Eqs.(62)-(64), (33) and (34). K here is 
the absolute minimum of %P- 
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Fig. 7. — The behaviour of Poynting jets about the equatorial boundary layer. The plots 
show: magnetic field structure (central frame, from Eq.(118)); the normalized density n e 
= ^j^- in the equatorial plane vs. the distance along the equatorial plane (upper right); 
the total proper current integrated about the equatorial plane vs. the distance along the 
equatorial plane (lower right); the ratio of the density to its equatorial value at the same r 
across the equatorial sheet (upper left) and the integrated proper current across the equa- 
torial sheet (lower left). The latter two quantities are plotted for some arbitrary value of r, 
r . The field structure (central frame) is plotted for -| = 0.9. The plots in the right panels 
are from Eqs.(106) and (107), while those in the left panels are from Eqs.(102), (103), (33) 
and (99) with v = \. 



